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We determine the third Zemach moment of hydrogen ({r 3 )(2)) using only the world data on 
elastic electron-proton scattering. This moment dominates the O(Za) hadronic correction to the 
Lamb shift in muonic atoms. The resulting moment, (r s }(2) = 2.71(13) fm 3 , is somewhat larger 
than previously inferred values based on models. The contribution of that moment to the muonic 
hydrogen 2S level is -0.0247(12) meV. 

PACS numbers: 36.10.Dr, 13.40.Em, 13.60.Fz 



I. INTRODUCTION 

Recent improvements in experimental techniques have 
led to the measurement of the 1S-2S interval in hydro- 
gen to an unprecedented accuracy of 2 parts in 10 14 . At 
the same time the QED corrections (expressed as pow- 
ers of a, the fine-structure constant) to this interval have 
been calculated through 0(a 7 ), with an important sub- 
set of 0(a 8 ) terms determined in the past few years (see 
Ref . for a comprehensive review of the theory and Ap- 
pendix A of for very recent developments). With 
these improvements the difference between experiment 
and QED theory (including recoil corrections) is domi- 
nated by hadronic size corrections, which affect the 1S-2S 
interval in the 10th significant figure. Uncalculated QED 
corrections probably affect the 12th significant figure Q. 

This situation creates both an opportunity and a prob- 
lem. The opportunity is to use the high precision mea- 
surements involving S-states to determine the rms charge 
radius of the proton at levels of precision of roughly 
1%. Recent values obtained in this way are: {r 2 }^ 2 — 
0.883(14) fmp, 0.891(18) fmQ, 0.869(12) frng, and 
0.875(7) fmQ. These consistent results reflect slightly 
different theoretical and experimental input. 

The murky history of experimental values for the pro- 
ton charge radius obtained from elastic electron-proton 
scattering data has recently been clarified by the compre- 
hensive analysis of That analysis of all of the world's 
data separated the charge from the magnetic scattering, 
incorporated (significant) Coulomb corrections 0, and 
carefully treated systematic (as well as random) uncer- 
tainties. The resulting value of (r 2 )p^ 2 — 0.895(18) fm is 
significantly higher than most older values, but is consis- 
tent with the atomic determinations. It is unlikely that 
new and relevant electron-scattering data will become 
available in the near future, and significant improvements 
over the 2% uncertainty are therefore unlikely during this 
period. 

The problem mentioned above is that the extremely 
precise measurements of hydrogen spectral lines cannot 



be used to test QED at anything approaching the level 
of accuracy of the QED theory, unless tricks are used 
to combine measurements of at least two experiments in 
such a way that the shortest-range hadronic and QED 
processes cancel between the measurements. Although 
this complementary approach is a very useful and ac- 
tive field of research 0, significant QED information has 
nevertheless been removed in the effort to eliminate the 
leading-order hadronic effects. Any alternative method 
that could provide a highly accurate value of the pro- 
ton charge radius would resolve much of this problem. A 
possible method is the measurement of the 2S-2P Lamb 
shift in muonic hydrogen, and such an experiment at PSI 
plans to determine the proton radius. 

The Lamb shift in electronic hydrogen is dominated 
by the (repulsive) radiative corrections on the electron 
line, which are much larger than the (attractive) vacuum 
polarization corrections on the photon line. The electron 
spends most of its time outside the polarization cloud in- 
duced in the electron Fermi sea. In the muonic-atom case 
the much smaller Bohr radius is within a significant por- 
tion of that cloud and the (electron) vacuum polarization 
dominates the QED corrections. The smaller radius also 
means that the hadronic size corrections are significantly 
more important, as well. The goal of the PSI experiment 

is to determine the proton charge radius, (r 2 )p^ 2 , to 1 
part per thousand. This requires a theoretical accuracy 
significantly better than .008 meV, which is the uncer- 
tainty introduced by a part per thousand error in the 
proton radius. The necessary theoretical developments 
have been recently reviewed in 01 1 Hi j and . 

Reference ^(| estimates that uncalculated QED dia- 
grams are likely to be smaller than 0.002 meV (although 
Ref. ^3 points out that the light-by-light-scattering con- 
tribution has not yet been calculated). The largest uncer- 
tainties are from hadronic contributions. Most of the lat- 
ter corrections fortunately are proportional (or nearly so) 
to the mean-square charge radius of the proton. All such 
terms can be simultaneously fit to the observed muonic 
Lamb shift together with the leading-order size correction 
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in Eqn. (1) below. One contribution, however, is signif- 
icantly different, and estimates of its size have shown 
considerable variation. We address this quantity in the 
next section. 



II. ZEMACH MOMENTS 

The primary hadronic size corrections are the Coulomb 
correction of 0(Za) A 1 the one-loop correction of O(Za) 5 , 
the two-loop correction of O(Za) 6 , and the hadronic-size 
modifications of the radiative and vacuum polarization 
corrections. In addition there are hadronic vacuum polar- 
ization corrections and polarizability corrections (which 
are significant, but which we ignore). 

The dominant and long- known |12j hadronic size cor- 
rection arises from modifying the Coulomb potential with 
the hadronic charge distribution, p^. Higher-order con- 
tributions of this mechanism can be obtained as well, and 
these were calculated many years ago in the context of 
muonic atoms 0, Q|. The first three orders of correc- 
tions for the nth S-state can be written in the general 
form: 

A£„ = |zh<M0)| 2 (V) - ^ 3 >( 2 ) 

+ {ZafF REL + {Zap) 2 F NR + •••), (1) 

where Z is the nuclear charge, (r m ) is the mth moment 
of the nuclear charge distribution (normalized to unit 
charge), p is the reduced mass of the muon- nucleus sys- 
tem, (j) n (0) is the muon wave function at the origin, and 
the Zemach moment Il5| (f 3 )(2) is defined by 

(r m )(2)= Jd 3 rr m p (2) (r), (2) 

where the convoluted (Zemach) charge density is given 
by 

P(i){r)=J d 3 z Pch (\z~r\) Pch (z) . (3) 

The nonrelativistic term -Fnr is part of the Coulomb 
correction of relative order (to the leading-order term 
in Eqn. (1)) (Za) 2 p 2 R 2 , where R is a generic proton 
radius, while the corresponding relativistic correction is 
determined by -Frel, and is of relative order (Za) 2 . 

The finite-size radiative corrections are discussed in de- 
tail in Refs. 0, E [ill and should scale like (r 2 ) p . The 
complete hadronic correction of 0(Za) e has never been 
worked out but the Coulomb approximation to it is small 
(see Refs. E3, El and above) and this is likely to be 
adequate. The remaining term is the one- loop contri- 
bution of O(Za) 5 , discussed in Refs. QUEl El E3 > 
which we discuss next. 

The one-loop correction of O(Za) 5 was considered by 
Pachucki|l6|. and was expressed in terms of the proton's 
Dirac form factors, F\ and F 2 . He also showed in the 



limit that the proton mass becomes very large that this 
contribution approached a relatively simple expression 
involving the same limit of the Sachs electric form factor, 
Ge(q 2 )- The expression 

(r 3 ) (2) = - r%G E (, 2 )-l + 9 V) p /3) (4) 
~ Jo 1 

can be verified by writing the form factor G^(q 2 ) as the 
Fourier transform of p c h(f) and repeatedly integrating 
by parts in Eqn. (4). Substituting Eqn. (4) for (^ 3 )(2) m 
Eqn. (1) verifies that Pachucki's simplified expression is 
just the O(Za) 5 term in Eqn. (1). The terms that vanish 
in this approximation can be considered as higher-order 
recoil corrections, as was noted in Refs. 0E3- According 
to the Coulomb approximation is good to within 10% 
for a simple model of the form factors. 

The convoluted density p( 2 ) (r) arises because each of 
the Coulomb photons is modified at short distances by 
the nuclear charge distribution. An alternative descrip- 
tion is that the modified Coulomb potential changes the 
muon wave function at short distances, and this has an 
effect on all expectation values^jl IT5L ITs[ . Thus the 
(?" 3 )(2) term above bears the same relationship to the 
full elastic one-loop contribution to the Lamb shift as 
the (traditional Zemach moment) (f)(2) term does to the 
full elastic one-loop contribution to the hyperfine struc- 
ture. We also note that analytic results exist for (r 3 )(2) 
for three simple charge distributions, including the dipole 
form factor case (viz., an exponential charge distribution, 
for which (r 3 ) (2) = 35V3(r 2 ) 3 / 2 /16)0. The necessity 
to resort to models of the proton form factor has lead to 
significant variations in the size expected for (r 3 )( 2 ). For 
a treatment using heavy-baryon effective field theory see 
Ref. [3. 

III. CALCULATIONS 

In Ref. the world data on electron-proton scatter- 
ing for momentum transfers q < 4 fm _1 have been ana- 
lyzed (for references to the data see @). The electric and 
magnetic Sachs form factors Ge(q 2 ) and Gm(<? 2 ) have 
been parameterized using a Continued Fraction (CF) 
expansion. It has been shown that this CF-expansion 
is more suitable than other parameterizations used in 
the past, and the contribution of the model dependence 
due to this choice has been evaluated. The longitudi- 
nal/transverse separation then is done during the global 
fit of the cross sections, an approach that is superior 
to the L/T-separations performed when determining Ge 
and Gm from individual data sets. 

The fit cross sections have been calculated from Ge 
and Gm including, in second-order Born approximation 
Q, the Coulomb distortion of the electron waves; this 
correction, although neglected in almost all analyses in 
the past, is important at low q. 

The data have been fitted using their random errors, 
and the error propagation treated via the error matrix. 
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The systematic uncertainties of the data have been taken 
into account by changing the data sets by the quoted er- 
ror, refitting and adding all resulting changes quadrat- 
ically, hereby obtaining a very conservative estimate of 
the systematic uncertainty of (r 3 )( 2 )- 

The result: (r 3 )(2) amounts to 2.71(13) fm 3 , where the 
error bar includes both random and systematic errors of 
the data, the latter dominating by far. For comparison 
we quote (r 3 )(2) for the "standard" dipole parameteriza- 
tion, which corresponds to an rms charge radius of 0.811 
fm and produces (r 3 )(2) =2.02 fm 3 . If the charge radius 
of the dipole model is scaled to 0.895 fm, the dipole result 
becomes 2.72 fm 3 , which is in very good agreement with 
the value we determined directly from the data. 

We should perhaps add a comment on the integral in 
Eqn. (4), which seems to indicate that, as a consequence 
of the l/g 4 -factor, the (r 3 }(2) depends on data at ex- 
tremely small q. It must be noted, however, that the 
Ge(0) term cancels the "1" at q ~ 0, and the q 2 (r 2 ) p /3 
term cancels the first, <7 2 -dependent term in a power se- 
ries expansion of GE(q)- Sensitivity studies have shown 
that the main contribution to the integral comes from 
the region q ~ 1.1 ± 0.5 fm -1 where the data base for 
electron-proton scattering is very good. 

We have also looked at the effect of two-photon ex- 
change (beyond Coulomb distortion) , which recently has 
been studied pflj in connection with differences between 
values of Ge extracted from Roscnbluth-scparations of 
cross sections and polarization transfer measurements. 
These corrections turn out to have a very minor effect; 



the value of (r 3 )(2) is increased by 0.02 fm 3 , i.e. by only 
a small fraction of the error bar. 



IV. CONCLUSIONS 

We have calculated the third Zemach moment for the 
charge distribution of hydrogen from the world's electron- 
proton scattering data. That moment is (r 3 )( 2 ) = 
2.71(13) fm 3 , which contributes -0.0247(12) mcV to the 
2S state of muonic hydrogen. For comparison the re- 
sult for the "standard" dipole model is an energy shift 
of —0.0185 meV. Our result is model independent, and 
removes the concerns of Refill]] > wn0 rioted that two 
simple models of the charge distribution with the same 
charge radius produced differences of 0.002 meV to the 2S 
energy shift. Our energy shift is somewhat larger than 
most recent values and significantly larger than a few. 
Our calculation removes most of the uncertainty from 
the contribution of the O(Za) 5 (elastic) finite-size term. 
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